Astronomy 8l Astrophysics manuscript no. 17206 © ESO 201 1 

July 15,2011 



^ 






Dissecting a hot molecular core: The case of G31.41+0.31'^ 

R. Cesaroni^ M.T. Beltran\ Q. Zhang^, H. Beuther^ and C. Fallscheer^'^ 

^ Osservatorio Astrofisico di Arcetri, INAF, Largo E. Fermi 5, 1-50125 Firenze, Italy e-mail: cesa@arcetri . astro . it 
^ Harvard-Smithsonian Center for Astrophysics, Cambridge, MA 02138, U.S.A. e-mail: qzhang@cf a. harvard, edu 
^ Max-Planck-Institut fiir Astronomie, Konigstuhl 17, D-69117 Heidelberg, Germany e-mail: beuther@mpia-hd.mpg . de 
^ Department of Physics and Astronomy, University of Victoria, 3800 Finnerty Road, Victoria, BC V8P 5C2, Canada e-mail: 
Cassandra . f allscheer@nrc-cnrc . gc . ca 

Received date; accepted date 

ABSTRACT 

o" 

(S| ■ Context. The role of disks in the formation of high-mass stars is still a matter of debate but the detection of circumstellar disks around 

0-type stars would have a profound impact on high-mass star formation theories. 

^ ' Aims. We made a detailed observational analysis of a well known hot molecular core lying in the high-mass star-forming region 

G3 1.41 +0.31. This core is believed to contain deeply embedded massive stars and presents a velocity gradient that has been interpreted 
either as rotation or as expansion, depending on the authors. Our aim was to shed light on this question and possibly prepare the ground 
for higher resolution ALMA observations which could directly detect circumstellar disks around the embedded massive stars. 
Methods. Observations at sub-arcsecond resolution were performed with the Submillimeter Array in methyl cyanide, a typical hot 
molecular core tracer, and ^^CO and ^^CO, well known outflow tracers. We also obtained sensitive continuum maps at 1.3 mm. 
Results. Our findings confirm the existence of a sharp velocity gradient across the core, but cannot confirm the existence of a 
bipolar outflow perpendicular to it. The improved angular resolution and sampling of the uv plane allow us to attain higher quality 
channel maps of the CH3CN lines with respect to previous studies and thus significantly improve our knowledge of the structure and 

r"| \ kinematics of the hot molecular core. 

O |. Conclusions. While no conclusive argument can rule out any of the two interpretations (rotation or expansion) proposed to explain 

I ■ the velocity gradient observed in the core, in our opinion the observational evidence collected so far indicates the rotating toroid as 

P ' the most likely scenario. The outflow hypothesis appears less plausible, because the dynamical time scale is too short compared to 

that needed to form species such as CH3CN, and the mass loss and momentum rates estimated from our measurements appear too 

^ . high. 

d . 
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> 

—1. ■ 1. Introduction from these OB-type stars are not sufficient to destroy the disk, 

QQ , which eventually turns into an ioni zed, rotating accretion flow 

^ High-mass stars are usually defined as those exceeding ~8 Mo, ^.^^^^ j^ j^e star (Sollins et al. l2005l Keto|2007l(. 

. based on the fact that stars above this mass limit do not have 
r^ ■ a pre-main sequence phase (Palla & Stabler [1993. This means Por these reasons it seems established that circumstellar ac- 
O ■ that accretion is ongoing until the star ignites hydrogen bum- cretion disks play a crucial role in the formation of all stars 
^^ ing and reaches the zero-age main sequence (ZAMS). At this and not only solar-type stars. This theoretical result contrasts 
Tl point the strong radiation pressure may halt and even reverse in- with the limited observational evidence of disks in high-mass 
> fall and thus stop further growth of the stellar mass. This led to (proto)stars. Only in recent years the number of disk candi- 
the so-called "radiation pressure problem". Recent studies have dates associated with luminous young stellar objects (YSOs) 
demonstrated that this limitation holds only in spherical symme- has significantly increased, mostly owing to the improvement 
h . try. As first envisioned by Nakano (fT987ll and recently demon- of (sub)millimeter interferometers in terms of angular resolution 
■ - strated by Krumholz et al. (2009) and Kuiper et al. (2010), ac- and sensitivity. The main problems that one has to face m this 
cretion through a circumstellar disk can explain the formation type of search are the large distances to the sources (typically 
of stars up to the upper limit of the initial mass function, by a few kpc) and the confusion caused by stellar crowding (OB- 
allowing part of the photons to escape along the disk axis and type stars form in clusters). These factors may explain the fail- 
boosting the ram pressure of the accreting gas through the small "^e to detect disks in 0-type stars, as opposed to the number of 
disk solid angle. It also appears that the powerful ionizing fluxes detections obtained for B-type stars (Cesaroni et al. 2007). In 

association with the most luminous YSOs one finds only huge 

~ Z ^ . , , , n ^ • -1 (^0.1 pc), massive (a few 100 Mq) cores, with velocity gradients 

Sena offprint requests to: R. Cesaroni, e-mail: ^. . —, > • ... •.„.•>. 

cesaOarcetri astro it suggesting rotation. These objects, named toroids , are likely 

* Based on observations carried out with the Submillimeter Array, non-equilibrium structures, because the ratio between the accre- 

The Submillimeter Array is a joint project between the Smithsonian Ho" time scale and the rotation period is shorter than for disks: 

Astrophysical Observatory and the Academia Sinica Institute of this implies that the toroid does not have enough time to adjust 

Astronomy and Astrophysics and is funded by the Smithsonian its structure to the new fresh material falling onto it (Cesaroni et 

Institution and the Academia Sinica. al.i2007i Beltran et al.i2011j). 
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With this in mind, one can see that understanding the forma- 
tion of the most massive stars (> 20 M©) may benefit from a de- 
tailed investigation of toroids, also because these objects might 
be hosting true circumstellar disks in their interiors. Moreover, 
their mere existence may set tighter constraints on theoretical 
models. Despite the number of candidates, the existence of ro- 
tating toroids is still a matter of debate. What is questioned is the 
nature of the velocity gradient, which is sometimes interpreted 
as expansion instead of rotation (see e.g. Gibb et al. 120041 and 
Araya et al. 2008 and references therein), thus suggesting that 
one might be seeing a compact bipolar outflow rather than a ro- 
tating core. In an attempt to distinguish between these two pos- 
sibilities and more in general to shed light on these intriguing 
objects, we have focused our attention on one of the best ex- 
amples: the hot molecular core (HMC) G31.41-h0.31 (hereafter 
G31.41). 



This prototypical HMC is located at a kinematic distance of 
7.9 kpc and was originally imaged in the high-excitation (4,4) in- 
version transition of ammonia (Cesaroni et al. 1994a) and in the 
(6-5) rotational transitions of methyl cyanide (CH3CN; Cesaroni 
et al. |1994bl) . The latter showed the existence of a striking veloc- 
ity gradient (centered at an LSR velocity of -96.5 km s"^) across 
the core in the NE-SW direction, already suggested by the dis- 
tribution and velocities of OH masers (Gaume & Mutel 1987). 
Follow-up interferometric observations with better angular res- 
olution and in high-energy tracers have confirmed this result 
and revealed the presence of deeply embedded YSOs, which in 
all likelihood explains the temperature increase toward the core 
center (Beltran et al. 2004, 2005, hereafter BEL04 and BEL05; 
Cesaroni et al. |2QTq1) . The G31.41 HMC is separated by ~5'' 
from an ultracompact (UC) Hii region, and overlaps in projection 
on a diff'use halo of free-free emission, possibly associated with 
the UC Hii region itself (see e.g. Fig. 2c of Cesaroni etaL [1998). 
The SPITZER/GLIMPSE images (Benjamin et al. l2003l) show 
that the HMC lies in a complex pc- scale region where both ex- 
tended emission and multiple stellar sources are detected. All 
these facts complicate the interpretation of the HMC structure 
and its relationship with the molecular surroundings and call for 
additional high-quality observations. 



An important aspect of a proper interpretation of the veloc- 
ity gradient observed in the HMC is the existence of an asso- 
ciated molecular outflow. A reasonable assumption, suggested 
by the analogy with low-mass YSOs, is that a disk actively un- 
dergoing accretion must be associated with a bipolar jet/outflow 
expanding along the disk rotation axis. Therefore, a way to dis- 
tinguish between the two interpretations of the HMC velocity 
gradient (expansion vs rotation) is to search for a molecular out- 
flow perpendicular to the gradient. Previous interferometric ob- 
servations in the ^^CO(l-O) line have suggested the presence of a 
collimated outflow directed SE-NW (Olmi et al. .1996} , but the 
poor uv coverage, especially on the shortest baselines, makes 
this result questionable. We have therefore conducted new inter- 
ferometric observations at 1.3 mm with the Submillimeter Array 
(SMA) in the typical HMC tracer CH3CN and in standard out- 
flow tracers such as ^^CO and ^^CO. In order to be sensitive to 
extended structures filtered out by the interferometer, we have 
also mapped the region over ~2' with the IRAM 30-m telescope. 
The observational details are given in Sect. [21 while the results 
are illustrated in Sect. [3] and discussed in Sect. IH Finally, the 
conclusions are drawn in Sect. [51 



Table 1. Parameters of SMA images. The RMS of the lines is 
estimated in each channel 
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2. Observations and data reduction 

Our observations have been performed with the SMA interfer- 
ometer and IRAM 30-m telescope. The former was used with the 
aim to image the ^^CO and ^^CO (2-1) lines and the CH3CN(12- 
11) transitions, but the 2 GHz bandwidth in both the LSB and 
USB allowed us to cover a much larger number of lines from a 
plethora of diff'erent molecules. Despite the richness of the spec- 
trum, in this article we will present only the results obtained for 
the above mentioned lines, which are the best tracers for our 
purposes. The 30-m telescope was used to fill the zero-spacing 
of the SMA maps in the ^^CO and ^^CO lines and thus help es- 
tablishing the morphology of the emission at all velocities and 
on all scales. Technical details of the observations are given in 
the next two sections. 



2.1. SMA Interferometer 

Obser vation s of G31.41 were carried out with the SMA [3 (Ho 
et al. 119961) in the 230 GHz band in the compact and ex- 
tended configurations. The correlator was configured to a uni- 
form spectral resolution of 0.41 MHz (0.6 kms"^) over the 
entire 4 GHz spectral window. With IF frequencies of 4 to 6 
GHz, the observations covered the rest frequencies from 219.3- 
221.3 GHz in the lower side band (LSB), and 229.3-231.3 
GHz in the upper side band (USB). The compact configura- 
tion data were obtained on 2007 July 09 with eight antennas. 
With 225 GHz zenith opacity of 0.05 to 0.1, the typical dou- 
ble side-band system temperatures were around 200 K. The very 
extended data were obtained on 2007 May 21 with seven an- 
tennas. At 225 GHz the zenith opacity was about 0.2 and the 
double side-band system temperatures were around 400 K. For 
both tracks, we used Vesta and 3C 273 for flux and bandpass 
calibrations. The time-dependent gains were calibrated using 
1751-h096 and 1830-h063. The phase center of the observations 
was a(J2000)=18M7^34^.315, ^(J2000)=-0r 12^45^9. The pri- 
mary beam of the 6-m antennas is about 55'' at the operating 
frequencies. 

The visibility data were calibrated in the IDL superset MIR0 
and MIRIAD (Sault et al. l2004l) . and were exported to MIRIAD 
format for imaging. The projected baselines of the combined vis- 
ibilities from both configurations range from 12 m to 500 m. The 
continuum data were constructed from line free channels and the 
continuum was subtracted from the line uv data. 

Continuum and channel maps of the CH3CN and CHs^^CN 
lines were created and cleaned in MIRIAD by weighting the data 



^ The Submillimeter Array is a joint project between the Smithsonian 
Astrophysical Observatory and the Academia Sinica Institute of 
Astronomy and Astrophysics, and is funded by the Smithsonian 
Institution and the Academia Sinica. 

^ The MIR cookbook can 

https ://w w w. cf a. harvard. edu/~ cqi/mircook. html 
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with "robust=0" to find a compromise between angular resolu- 
tion and sensitivity to extended structures. The resulting synthe- 
sized beam is approximately (X/ 88 x (7/75, with P. A. 53°. The 1 cr 
rms is 5.8 mJy/beam in the continuum, and 89 mJy/beam in the 
line images with 0.6 km s"^ spectral resolution. The ^^CO and 
^^CO (2-1) line data were imaged in a similar fashion, but using 
natural weighting to enhance extended emission. The resulting 
synthesized beam is a.'89 x 1773 with P.A. 53°. The 1 cr rms in 
the 1 km s"^ channel maps is 50 mJy/beam. 

The ^^CO and ^^CO data were also combined with the 30- 
m data (see Sect. 12.21) to recover extended structure resolved 
out by the interferometer. The single-dish maps were Fourier- 
transformed and then suitably sampled in the uv domain, then 
merged with the SMA data. Finally natural weighted maps were 
created using the same procedure as for the non-merged data. 
The resulting synthesized beam is 2! 3 x V!3 with P.A. 63° and 
the 1 cr rms in the 1 km s"^ channel maps is 66 mJy/beam. 

2.2. IRAM 30-m telescope 

Maps of G31.41 with the 30-m telescope were made on 
November 4th, 2007. The HERA multi-beam receiver was used 
to cover a region of 2' x2' centered on »( J2000)= 1 8M7"'34'.3 and 
^(J2000)=-01°12'45r9, using the on-the-fly-mode. To prevent 
systematic eff'ects, the region was scanned alternatively along 
right ascension and declination, and all the data were eventu- 
ally averaged to obtain the final map. The receiver was tuned 
to the frequencies of the ^^00(2-1) (HERAl) and ^^C0(2- 
1) (HERA2) lines, which were hence observed simultane- 
ously. Both lines were covered with the VESPA autocorrela- 
tor, with 0.4 km s"^ spectral resolution. The reference position 
used for all maps was a(J2000)= 18^5^15^.24 and ^(J2000)=- 
00°55'56V09. This was carefully chosen from CO surveys such 
as that by Sanders et al. ( 1986), and checked to be free of emis- 
sion in the observed lines. Data were reduced with the pro- 
gram CLASS of the GILDAS packageEl and channel maps with 
1 km s"^ resolution were created. 



3. Results 

In the following we illustrate the results obtained, first for the 
typical hot core tracers (continuum, CH3CN, and CHs^^CN 
linefl ) and then for the outflow tracers (^^CO and isotopomers). 
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Fig. 1. Spectra obtained by integrating the line emission over the 
whole core. Top. Ground-state ^-components of the /=12^11 
transition of CH3CN (numbers above the spectrum) and of the 
CHs^^CN isotopomer (numbers below the spectrum). Note how 
most /sT-lines are blended with each other and/or with lines of 
other molecules. Middle. Vibrationally excited (v8=l) lines of 
CH3CN. The KJ quantum numbers of the transitions are indi- 
cated above the spectrum. Bottom. Ground state K=3 compo- 
nent. The dotted lines mark the velocity ranges used to produce 
the maps in Figs. [2] and O 



3.1. Core tracers 

Compared to previous Plateau de Bure interferometer (PdBI) ob- 
servations in the same lines (BEL04, BEL05), our SMA observa- 
tions benefit from a better uv coverage and more circular beam, 
thus allowing us to attain sub-arcsec resolution in all directions 
(see Table [H for comparison, the synthesized beam of the PdBI 
images was irixa.'5 with PA=189°). 

To give an idea of the complexity of the line emission from 
the G31.41 HMC, in Fig. [T] we show a spectrum covering all 
CH3CN(12-1 1) transitions with ^ < 10, obtained by integrating 
the emission over the whole core. For an estimate of the CH3CN 
and CHs^^CN line parameters we refer to the study of BEL05 
(see their Table 4). We did not attempt the identification of all 
lines detected, because this goes beyond the purposes of the 



available 



^ The GILDAS software h 

|http://www.iram.fr/IRAMFR/GILDAS 

^ Parameters of the CH3CN and CHs^^CN transitions can be found, 
e.g., in the "splatalogue" database: http://www. splatalogue.net/ 1 



at 



present study. However, we note that a number of CH3CN and 
CHs^^CN transitions are heavily blended with each other and/or 
with other lines. Only the CH3CN ^=2,3,4,8 and the CHs^^CN 
K=2,6 components appear free enough of contamination to be 
used for a detailed investigation. A similar consideration holds 
for the vibrationally excited (v8=l) lines, shown in the same fig- 
ure. Here, only the (^=3,/=-l) and (^=6, /=1) transitions are 
clearly detected and sufficiently separated from other lines to 
be considered for our study. Note that our bandwidth covers all 
V8=l lines down to the (11,-1) at 220788.016 MHz, but the lines 
below 221 . 19 GHz are too weak or too blended: this is the reason 
why the CH3CN( 12-11) V8=l spectrum in Fig. [T] is displayed 
only above this frequency. 

Figure [21 shows continuum and line emission maps obtained 
by averaging the emission under the line and Table [2] lists the 
main parameters of the continuum emission. The most obvi- 
ous feature is that the maps of all tracers are slightly elongated 
in the NE-SW direction, suggesting that the distribution of the 
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Fig. 2. Maps of the line and continuum emission toward the G31.41 HMC. The line maps have been obtained by averaging the 
emission under the whole line profile, namely from 88.9 km s"^ to 103.9 km s"^ (see Fig.[T] bottom). The tracer used is indicated 
in the top right of each panel. The two white points mark the positions of the free-free-continuum sources detected by Cesaroni 
et al. (|2010). The contour levels start from 0.072 Jy/beam and increase in steps of 0.12 Jy/beam for all line maps, and range from 
0.035 Jy/beam to 1.785 Jy/beam in steps of 0.175 Jy/beam for the continuum map. Off'sets are measured with respect to the phase 
center a(J2000)=18M7"^34'.315 ^(J2000)=-0r 12^45^9. The ellipse in the bottom right denotes the FWHP of the synthesized beam. 



molecular gas is flattened approximately in the direction defined 
by the two free-free continuum sources detected by Cesaroni et 
al.(l20T0l). 

From the same figure one can also note that the high-energy 
transitions and the isotopomer emission peak approximately at 
the same position as the continuum, whereas the maximum of 
the low-energy K=2 component appears to be slightly off'set 
from it. Such an off'set is evident for all components with K < 6, 
while the K=l and 8 line maps look more circular and peak at 



the core center. Nothing can be said about the K=9 transition, 
which is heavily blended with the ^^C0(2-l) line. A plausible 
interpretation is that the opacity decreases with increasing exci- 
tation energy. In confirmation of that, one can see that this eff'ect 
is even more prominent in the maps obtained by integrating the 
line emission over a narrow velocity range around the systemic 
LSR velocity, i.e. from 95.5 to 97.3 km s"^(see Fig. [3]). Now, the 
K=2 emission clearly splits into two peaks located on opposite 
sides with respect to the continuum and also the corresponding 
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Fig. 3. Same as Fig.[2l with the line maps obtained by averaging the emission from 95.5 km s"^ to 97.3 km s"^ (see Fig.[TJ bottom). 
Contour levels start from 0.15 Jy/beam and increase in steps of 0.2 Jy/beam for all line maps, and range from 0.035 Jy/beam to 
1.785 Jy/beam in steps of 0.175 Jy/beam for the continuum map. 



Table 2. Continuum parameters, i.e. peak position, peak bright- 
ness temperature in the synthesized beam, integrated flux den- 
sity, full width at half maximum (FWHM), and deconvolved an- 
gular diameter 
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CHs^^CN line map shows a similar pattern. We conclude that 



opacity plays a crucial role in the CH3CN emission of this HMC 
and cannot be neglected in the analysis of the line emission. 

A qualitative impression of the CH3CN emission can be ob- 
tained from visual inspection of the channel maps shown in 
Fig. m In order to improve the signal-to-noise, such maps have 
been obtained by averaging the emission in the K=2, 3, and 4 
components. These transitions have been chosen because they 
do not seem to be afl'ected by significant contamination by other 
lines and have comparable excitation energies, which suggests 
that they are likely tracing the same gas. 

Clearly, the emission gradually drifts from SW to NE with 
increasing velocity and outlines an "8-shaped" pattern close to 
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Fig. 4. Channel maps of the CH3CN(12-11) emission obtained by averaging the K=2, 3, and 4 components. The number in the 
top left of each panel gives the corresponding LSR velocity. The values of the contour levels are indicated by the tick marks in the 
intensity scale to the right. The cross denotes the phase center. The ellipse in the bottom right of the last panel is the full-width-at-half 
power of the synthesized beam. 



the systemic velocity of ~96.5 km s"^. These characteristics re- 
semble those observed in rotating circumstellar disks around 
low-mass stars (see e.g. Fig. 1 of Simon et al. 2001 ). We will 
discuss this problem in Sect. 14.31 A more quantitative analysis 
of the velocity field traced by the CH3CN emission is attained 
- in analogy with the previous study by BEL04 - by fitting the 
^=0 to 4 components simultaneously, after fixing their separa- 
tions in frequency to the laboratory values and forcing the line 
widths to be equaQ. This procedure improves the quality of the 
fit with respect to fitting a single K line, but we stress that a simi- 
lar result is obtained by fitting a single Gaussian to any of the un- 
blended K components. A map of the LSR velocity is presented 
in Fig. [21 which confirms the existence of a clear velocity gra- 
dient directed approximately NE-SW (RA.^68°) and roughly 
centered on the peak of the continuum emission. A diff'erent and 
more detailed representation of the same velocity gradient is pre- 
sented in Fig.O where the position-velocity plots along the ve- 
locity gradient (P.A.=68°) and perpendicular to it (P.A.=-22°) 
are shown for three diff'erent lines of CH3CN and CHs^^CN. 
Note that to enhance the signal-to-noise ratio in the plots, we 
have averaged the emission along the direction perpendicular to 
the cut along which the position-velocity plot is calculated. As 
expected, no velocity trend is seen in the plots along P. A. =-22°. 



It is worth comparing our results to those of BEL04. With re- 
spect to their Fig. 2e, our Fig.[5]presents a slightly broader veloc- 
ity range (~5.5 km s"^ instead of ~4 km s"^ ), over a larger region 
(~3'' instead of -2''). These diff'erences are probably caused 
by the diff'erent angular resolutions and sensitivities. A more 
evident discrepancy appears in the map of the CHsCNC 12-11) 
emission, i.e. our Fig. [3] and their Fig. 1. In the former, the K=2 
emission appears to outline two peaks located approximately to 
the NW and SE with respect to the core center, whereas in Fig. 1 
of BEL04 the peaks lie to the E and W. However, one should 
keep in mind that this map was obtained integrating the emis- 
sion under the ^=0,1, and 2 components. Moreover, we have 
averaged the emission over a narrow interval around the peak, 
whereas BEL04 integrated under the whole line profile. Finally, 
the angular resolution is diff'erent for the two data sets. In order 
to attain a more consistent comparison, we have reconstructed 
our maps using the same clean beam as BEL04 and integrating 
over the same frequency interval. The new map (shown in Fig. [7]) 
is now much more similar to that of BEL04, demonstrating how 
the apparent morphology of the emission may depend signifi- 
cantly on the resolution. Given the larger number of antennas 
and more circular beam, we believe that our SMA images repro- 
duce the structure of the emission more faithfully than the old 
maps by BEL04 and BEL05. 



^ For this purpose, program CLASS of the GILDAS package was 
used, after extracting the relevant spectra from the data cube and con- 
verting them to the appropriate format. 
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Fig. 6. Position-velocity plots of along P. A. =68° (left panels) and -22° (right panels) for three different lines (indicated in each 
panel). Offsets are measured with respect to the phase center. Contour levels start from 0.12 Jy/beam and increase in steps of 
0.12 Jy/beam. The cross in the bottom left of each box denotes the angular and spectral resolutions. 



3.2. Outflow tracers 

One of the purposes of the present study was to compare the 
structure and kinematics of the HMC with those of a possible 
bipolar outflow associated with it. The existence of such an out- 
flow had been suggested by Olmi et al. J1996b , whose Plateau 
de Bure images in the ^^CO(l-O) line seem to reveal two nar- 
row lobes oriented SE-NW (see their Fig. 5). With this in mind, 
we have imaged G31.41 in the ^^CO and ^^CO /=2-l rota- 
tional transitions with both the SMA and 30-m IRAM telescope, 
to recover also the emission filtered out by the interferometer. 
Figure [8] shows channel maps in the ^^C0(2-l) line and effec- 
tively demonstrates the complexity of the region. Only at high 



velocities the emission appears quite compact, whereas close to 
the systemic velocity the structure presents a complicated pat- 
tern. In particular, from 95 to 100 km s"^ the emission is uni- 
formly distributed over an extended region: this means that the 
line opacity is sufficiently high that one can see only the surface 
of the cloud at these velocities, consistent with the relatively low 
values of the brightness temperature (10-15 K). 

One arrives at the same conclusion by looking at Fig. [9l 
where a comparison is shown of the single-dish spectra obtained 
by averaging the emission of three different CO isotopomers 
over a square region, 30'' in size, centered at the HMC position. 
While the ^^CO and ^^CO profiles are double-peaked with a dip 
at -97 km s"^ the optically thin(ner) C^'^0(2-l) line (Cesaroni, 
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Fig. 8. Channel maps in the ^^C0(2-l) line obtained by combining single-dish (IRAM 30-m) and interferometric (SMA) data. The 
number in the top left of each panel gives the corresponding LSR velocity. The values of the contour levels are indicated by the 
tick marks in the intensity scale to the right. The ellipse in the bottom right of the last panel is the full-width-at-half power of the 
resulting instrumental beam. 
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Fig. 5. Overlay of the map of the 1.3 mm continuum emis- 
sion (contours) on that of the CH3CN(12-1 1) line velocity (color 
scale). Contour levels range from 0. 1 to 2 in steps of 0.2 Jy/beam. 
Offsets are measured with respect to the phase center. The two 
white dots denote the free-free continuum sources detected by 
Cesaroni et al. (120 101) . The synthesized beam is shown in the 
bottom left. 



unpublished data) is Gaussian and peaks right at the velocity 
of the dip in the other two lines. This is an indication of self- 
absorption and thus of high optical depth in the ^^CO and ^^CO 
transitions. 

An apparent feature of the ^^C0(2-l) profile is the presence 
of broad wings. From Fig. [8] one sees that this high- velocity 
emission originates from 2-3 compact structures located roughly 
to the E (red-shifted) and W (blue- shifted) of the HMC. Whether 
this morphology is consistent with that of a bipolar outflow can- 
not be trivially decided on the basis of the evidence presented 
so far and requires a detailed analysis and discussion that we 
postpone to Sect. 14.21 



4. Discussion 

The main purpose of the present study is to shed light on the na- 
ture of the velocity gradient observed in this HMC (e.g. Fig. [5]). 
As already mentioned, BEL04 and BEL05 favored the rotat- 
ing toroid scenario. In contrast, other authors (Gibb et al. 120041 
Araya et al. 120081) preferred the outflow interpretation, thus pos- 
ing a problem that we wish to address with our new SMA data. 
Besides the (possible) existence of rotation and/or outflow, 
the situation in G31.41 is complicated by the existence of infall, 
detected by Girart et al. (2009) as an inverse P-Cygni profile of 
the C^'^S(7-6) line observed with the SMA. Therefore, before 
making a comparative discussion of the diff'erent models for the 
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Fig. 7. Map of the CH3CN(12-1 1) emission obtained by aver- 
aging the emission under the ^=0 and 1 components. The image 
has been reconstructed with a clean beam (ellipse in the bottom 
right) equal to that of the images of BEL04, i.e. V!lxQf!5 with 
PA=189°. The contour levels range from 75 to 600 in steps of 
75 mJy/beam. Offsets are measured with respect to the phase 
center. 

velocity gradient, we need to verify to what extent our CH3CN 
measurements may be affected by infall. 

4.1. Evidence for infall 

As already discussed in Sect. 13.21 Fig. [9] clearly reveals that 
the ^^CO and ^-^CO (2-1) lines are affected by self-absorption. 
However, only a weak asymmetry is seen in the line profile, with 
the red- shifted peak being slightly more prominent than the blue- 
shifted one. If such an asymmetry were caused by inf ailing gas, 
the blue-shifted peak should be stronger; hence we can conclude 
that no obvious evidence of infall is present on the large scale 
(30'' or 1.1 pc) over which the spectra in Fig.[9]have been calcu- 
lated. It thus seems that the infall detected by Girart et al. (2009]) 
occurs close to the HMC. Do we see any evidence of this in 
the HMC tracers that we observed? All CH3CN and CHs^^CN 
(12-11) lines shown in Figure[T]have Gaussian profiles with no 
hint of (self-)absorption, unlike the C^^S(7-6) line observed by 
Girart and collaborators. However, this is not surprising, because 
both the excitation energy and critical density of the latter (49 K 
and ~ 10^ cm"^) are less than those of the former transitions 
(58-926 K and > 10^ cm'^). This implies that the C^^S(l-6) 
line traces large radii of the HMC that are characterized by a 
relatively low temperature, and thus absorb the hot continuum 
photons emitted from the central region of the core. This effect 
does not apply to the CH3CN(12-11) transitions, because these 
arise from a smaller hot shell, whose excitation temperature is 
much more similar to the brightness temperature of the dust con- 
tinuum. 

Despite the apparent lack of evidence for infall in the 
CH3CN transitions, we have attempted a more detailed analysis 
to see whether some hint of asymmetry was present in the line 
profiles. The idea is that if the line shape is skewed toward the 
blue (although very weakly), a Gaussian fit should give a peak 




100 
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Fig. 9. Spectra of the ^2C0( (solid histogram), ^^CO (dashed), 
and C^^O (dotted; Cesaroni, unpublished data) (2-1) lines ob- 
tained by averaging the emission over a square region 30'' in 
size, centered on the HMC. The data have been taken with the 
IRAM 30-m telescope. 



velocity less than the systemic velocity. With this in mind, in 
Fig. [To] we plot the line first moment (solid curves) as a function 
of distance from the HMC center along an axis with P. A. =-22°, 
i.e. perpendicular to the direction of the velocity gradient. For the 
sake of completeness, also the line zero moment (integrated in- 
tensity) along the same direction is shown (dashed curves). This 
is done for the ^=2,3,4, and 8 components of the ground- state 
transition and the (3,-1) line of CH3CN V8 = l. Evidently, up to 
^=8 the line velocity attains a minimum value close to the HMC 
center, whereas this minimum is not seen in the V8 = l line. This 
proves that the lower energy lines are slightly skewed to the blue 
toward the bright continuum peak, as expected for weak red- 
shifted absorption caused by infall. This effect is not seen in the 
highest energy line, coming from the innermost, hottest layers of 
the core. 

We conclude that some of the CH3CN( 12-11) lines are af- 
fected by red-shifted absorption, although in a much less promi- 
nent way than the C^^S(7-6) line measured by Girart et al. 
([2009). Note that this absorption causes only a marginal shift 
in the line velocity, of -0.5 km s"^ and is going to have only 
negligible effects on our study of the velocity shift of ~5 km s"^ 
observed across the HMC. 



4.2. The NE-SW velocity gradient 

The existence of a velocity gradient in the CH3CN lines is 
proved beyond any doubt by Figs. IH [3 and[6j The high angu- 
lar resolution allows us to establish that the velocity is shifting 
smoothly across the core, and one can reasonably exclude that 
this shift is caused by two unresolved velocity components (sub- 
cores) with different velocities. The question we intend to ad- 
dress here is whether this velocity gradient is caused by rotation 
of the HMC about a SE-NW axis, as hypothesized by BEL04, 
or to expansion in a bipolar outflow oriented NE-SW, as pro- 
posed by Gibb et al. (2004) and Araya et al. (2008). How can 
one distinguish between these two scenarios? As explained in 
Sect. [H it is common belief that disks and outflows are tightly 
associated, with the latter being ejected along the rotation axis 
of the former. Therefore, if the velocity gradient in the G31.41 
HMC is caused by rotation, one would expect to detect a bipolar 
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Fig. 10. Solid curves represent the first moments of sev- 
eral CHsCNC 12-11) components (indicated in the top left of 
each panel) estimated along a cut through the continuum peak 
(a(J2000)=18M7"^34^.31 ^(J2000)=-0r 12^4670) with RA .=- 
22°, i.e. perpendicular to the velocity gradient in Fig. [5] The 
dashed curves are the zero moments of the corresponding lines. 



outflow perpendicular to it, on a larger scale. In contrast, if the 
velocity gradient is tracing the "root" of a bipolar outflow, on 
a larger scale this outflow should become clearly visible along 
the same direction defined by the (small-scale) velocity gradient. 
Our combined single dish and interferometric observations in 
typical outflow tracers such as ^^CO and ^^CO should hence be 
well suited for our purposes, given the sensitivity to both small 
(a few arcsec) and large (a couple of arcmin) scales. 

The results reported in Sect. 13.21 are quite ambiguous. The 
blue- and red- shifted ^^CO and ^^CO gas is oriented like the 
CH3CN velocity gradient, suggesting that the two are manifes- 




tations of the same phenomenon. The existence of a NE-SW 
bipolarity in the ^^CO emission can be appreciated in Fig. [TTl 
where we show the blue- and red- shifted ^^C0(2-l) emission 
in pairs of channels equally ofl'set from the systemic velocity. 
One is tempted to conclude that ^^CO is indeed associated with 
a bipolar outflow, whose densest component is traced by the 
CH3CN emission. However, this bipolarity in the CO maps is 
seen only on a small scale, because no evidence of high- velocity 
CO emission is found beyond ~5'' (i.e. 0.2 pc) from the HMC, 
implying an unusually small size for a typical outflow associated 
with a high-mass star-forming region. This casts some doubt on 
the outflow interpretation. Below we discuss the two hypotheses 
(outflow and toroid) in detail. 



4.2.1. Bipolar outflow 

The outflow hypothesis is supported by a couple of facts. First 
of all, if the CH3CN gradient is caused by a rotating toroid, one 
expects to find a bipolar outflow perpendicular to it, but this is 
not seen in our maps. Second, the position velocity plot of the 
^^CO emission in the direction of the CH3CN velocity gradi- 
ent suggests that the gas velocity is proportional to the distance 
from the star (as discussed later in Sect. 14.2.21) , consistent with 
the Hubble-law expansion observed in molecular outflows from 
YSOs. 

In order to check the plausibility of the outflow scenario, 
we have calculated the outflow parameters from the CH3CN, 
CH3^^CN, ^^CO, and ^^CO lines, in the latter two cases using 
the combined 30-m and SMA data. These are given in Tabled 
where we report the mass of the outflow, M, the momentum, P, 
the energy, E, and the corresponding rates obtained by dividing 
the previous quantities by the dynamical time scale (see below). 
In practice, M = Y^iMi, P = Y^imiVi, and E = 2/(l/2)m/y^, 
where the sums are extended over only those channels, /, falling 
in the velocity ranges given in the footnotes of Table [3l and m/ 
is the mass moving with velocity Vt relative to the systemic ve- 
locity. This mass is computed by integrating the line emission 
inside the regions corresponding to the 5cr level of the blue- 
and red-shifted emission. Note that the velocity intervals were 
chosen by inspecting the channel maps and selecting only those 
channels where the emission was sufficiently strong and at the 
same time not aff'ected by the missing flux problem close to the 
systemic velocity. In our calculations we assume a temperature 
of -100 K for both molecules, intermediate between the peak 
brightness temperature (77 K) of the ^^C0(2-l) line and the 
rotational temperature (164 K) estimated by BEL05 from the 
CH3^^CN(6-5), CH3^^CN(12-11), and CH3CN(6-5) V8=l lines 
(see their Fig. 6). While this temperature may seem too high for 
the outflow component traced by the CO isotopomers, assuming 
50 K instead of 100 K would reduce the outflow parameters by 
only a factor 0.6. The abundance of CH3CN and ^^CO relative to 
H2 are assume d equa l to 10"^ and 10""^ respectively (see e.g. Van 
Dishoek et al. 119931) , while the isotopic ratios ^^CO/^^CO and 
CH3CN/CH3^^CN are taken equal to 50 after Wilson & Rood 
( 1994), for a galactocentric distance of 4.5 kpc. The dynamical 
time scale of the outflow, ^dyn - 4 x 10^ yr, is calculated from 
the maximum size (-0.12 pc) and velocity (~30 km s"^) of the 
^^CO lobes. Note that the (unknown) inclination of the outflow 
with respect to the line of sight has not been taken into account 
in our estimates. 

From Table [3] one notes that the values increase with de- 
creasing abundance of the molecule. This can be explained in 
terms of decreasing optical depth, because higher opacities lead 
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Fig. 11. Channel maps of the ^^C0(2-l) Hne emission observed with the SMA. Each box contains pairs of maps corresponding 
to the blue- (solid contours) and red-shifted (dashed) emission at the same velocity offset (in absolute value) from the the systemic 
velocity (-96.5 km s"^). The corresponding LSR velocities are indicated in the top right of the box. Contour levels increase from 
0.5 Jy/beam in steps of 0.5 Jy/beam. Also shown is a map of the CH3CN( 12-11) K-A line emission averaged over the range 
95.2-96.7 km s"^ . The ellipse in the bottom left denotes the synthesized beam of the observations. 



to an underestimate of the column density and hence of the mass. 
Indeed, from the ratio between the ^^CO and ^^CO emission in 
the line wings one derives opacities as high as -30 for ^^CO. 
Therefore, the most reliable estimates should be those obtained 
from CH3 ^^CN. These depend significantly on the CH3CN abun- 
dance, which is known to present considerable variations in 
molecular clouds. However, the value assumed in our calcula- 
tions is one of the highest found in the literature and we there- 



fore believe that our estimates are likely to be lower limits. This 
leads us to another consideration. Our estimates are very close 
to the largest outflow parameters ever measured (see e.g. Lopez- 
Sepulcre et al. 2009 and Wu et al. i2QQ4b . corresponding to a YSO 
powering the outflow of at least ~ 10^ L©. More precisely, us- 
ing Wu et al. (12004b relationship between bolometric luminosity 
(Lboi) of the powering source and outflow momentum rate, from 
the value in Table [3] one obtains Lboi = 6 x 10^ L©. This is much 
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Table 3. Outflow parameters calculated from diff'erent lines. The 
dynamical time scale is ^dyn - 4 x 10^ yr and a gas temperature 
of 100 K is assumed in the calculations 



parameter 


llQQa 


i^CO^ 


CH3CN 




M(Mo) 


3.8 


20 


60 


290 


P(Mokms-i) 


48 


230 


350 


1200 


£(Loyr) 


5.9 X 10^ 


2.2 X 10^ 


1.7 xlO^ 


4.2 X 10^ 


M(Moyr-i) 


9.5 X 10-4 


5 X 10-3 


0.015 


0.07 


P(Mokms-Vr"0 


0.012 


0.057 


0.087 


0.30 


E{L^) 


14.7 


54.6 


41.8 


104 



^ blue wing 
127.9 km s- 
^ blue wing 
112.9 km s- 
^ blue wing 
103.9 km s" 
^ blue wing 
101.5 km s" 



from 78.9 to 89.9 km s-\ red wing from 105.9 to 

^ Abundance relative to H2: 10-"^ 

from 77.9 to 88.9 km s-\ red wing from 103.9 to 

^ Abundance relative to H2: 2 x 10-^ 

from 88.9 to 91.9 km s-\ red wing from 100.9 to 

^ Abundance relative to H2: 10-^ 

from 90.1 to 92.5 km s-\ red wing from 99.1 to 

\ Abundance relative to Ho: 2 x 10-^^ 



higher than the luminosity estimated for the HMC (~ 10^ L©) 
by Osorio et al. C2009') and that obtained from the corresponding 
IRAS fluxes (2.6 x 10^ L©; see Cesaroni et al. 1994a). Moreover, 
the momentum rate of 0.3 M© km s"^ yr"^ is an order of magni- 
tude higher than that estimated by Cesaroni et al. (2010) from 
VLA observations of the free-free continuum emission (~0.03- 
0.06M©kms-iyr-^). 

It is also worth noting that the velocity gradient in G31.41 
seems to involve the whole core and not only the gas emitting 
in the line wings. Usually, in outflow sources the emission close 
to the systemic velocity - i.e. that inside the FWHM of the line 
- traces the molecular core; in contrast, in G31.41 most of the 
CH3CN emission is afl'ected by the velocity gradient, as one can 
see from the channel maps in Fig. |4l where only emission over 
a limited velocity range (~2-3 km s"^) appears to arise from the 
central region. As discussed below in Sect. 14.31 this situation is 
reminiscent of that observed in circumstellar disks around lower- 
mass YSOs. Our estimates of the outflow parameters based only 
on the line-wing emission are likely lower limits, which makes 
the case of G31.41 even more extreme compared to typical out- 
flows. 

In conclusion, the outflow hypothesis seems to yield values 
for the outflow parameters that are too high. In addition, one 
should note that the typical parameters given in studies such 
as those quoted above refer to single-dish observations of pc- 
scale flows, an order of magnitude higher than that observed in 
G31.41. Although we cannot exclude that in G31.41 one is ob- 
serving the earliest stages of the expansion, it is questionable 
that the parameters of a young, compact outflow are greater than 
those typical of much older, extended outflows. 

Finally, we note that the dynamical time scale of the putative 
outflow is an order of magnitude shorter than the time needed to 
form typical hot core species - such as methyl cyanide - accord- 
ing to theoretical models (see e.g. Charnley et al. 1992). Our 
estimate of ^dyn is afl'ected by large uncertainties, owing to the 
unknown inclination angle and the difl&culty in tracing the whole 
extent of the lobes. However, to increase ^dyn by a factor ^10 one 
has to assume that the outflow lies very close (^6°) to the line of 
sight and/or that we detect only the most compact part of lobes 
extending over a region 10 times larger than that imaged in our 
CH3CN and CO maps. The latter explanation is ruled out by our 



SMA-h30m combined maps that do not reveal any large-scale 
bipolar outflow. The former would imply a significant overlap 
between the blue- and red- shifted lobes in the plane of the sky, 
which is not seen in Fig. [141 Indeed, from the observed sepa- 
ration between the peaks of the red- and blue-shifted emission 
(~rO, assuming that the intrinsic length of the CH3CN lobes 
is comparable to the radius of the CH3CN core measured in the 
plane of the sky (-T'). one can obtain a rough estimate of the in- 
clination angle of arc sin(0. 5) = 30°. Correspondingly, ^dyn would 
increase by only a factor 1.7. 

Given the number of problems encountered in the outflow 
scenario, we investigate also a second possibility, namely that 
the velocity gradient of the HMC is due to a rotating toroid. 



4.2.2. Rotating toroid 

In an attempt to shed light on the velocity field in the HMC and 
its molecular surroundings, we overlay in Fig. [121 the position- 
velocity plots along the velocity gradient in high- and low- 
density tracers. For the sake of completeness we show both the 
SMA and the combined 30m-hSMA data in the ^^CO and ^^CO 
(2-1) transitions. Note that the cut was made along P. A. =68° and 
all plots were obtained after averaging the emission along the di- 
rection perpendicular to the cut. This has the twofold purpose of 
increasing the S/N of the plots and taking into account emission 
at all off'sets along the (putative) rotation axis. 

As already explained, the ^^CO and ^^CO emission close 
to the systemic velocity is highly opaque and does not convey 
any information on the HMC, but is instructive to study the 
high- velocity and/or low-density gas around it. Vice versa, the 
CH3CN emission is an excellent HMC tracer, but is not detected 
at high velocities and/or beyond ~2'' from the center. Combining 
all these tracers is the only way to perform a detailed and com- 
plete analysis of the velocity field in this core. Indeed, Fig. [121 
is very instructive. It shows that the velocity trend observed in 
CH3CN is complementary to that seen in the CO isotopomers 
and fills the gap caused by the interferometer resolving out the 
extended, bulk CO emission. 

A more thorough inspection of this figure reveals that, de- 
spite the similarity, the ^^CO and ^^CO plots are significantly 
diff'erent. While the ^^CO and CH3CN patterns match very well 
in the overlapping regions (i.e. between 89 and 94 km s"^ and 
between 102 and 105 kms"^), the ^^CO emission is off'set by 
~r' from the CH3CN and ^^CO patterns, to the SW in the blue- 
shifted part, and to the NE in the red- shifted part. An interest- 
ing diff'erence between CH3CN and ^^CO, on the one hand, and 
^^CO, on the other, is that the high-velocity emission is found 
at relatively large off'sets in ^^CO, whereas in ^^CO it appears to 
peak close to the HMC center. 

Is it possible to find a coherent interpretation of all these facts 
in the rotating toroid scenario? A possible explanation is that of 
a self-gravitating structure with a heavy, compact stellar cluster 
at the center. In the outer regions (traced by ^^CO), the gravita- 
tional field is determined by the gas mass and the rotation curve 
flattens, whereas close to the center (where ^^CO and CH3CN 
are detected) the velocity field resembles Keplerian rotation, be- 
cause the gas mass becomes comparable to the stellar mass. This 
scenario is analogous to that studied by Berlin & Lodato (1999), 
who determined the rotation curve of a self-gravitating disk with 
a central star. 

In order to check the viability of this interpretation, below 
we first estimate the mass of the HMC and then verify whether 
this suffices to support rotational equilibrium. 
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Fig. 12. Position-velocity plots along the direction with P. A. =68° passing through the HMC. The offset is measured from the 
phase center, positive toward NE. Contours correspond the CHsCNC 12-11) K=4 component emission, while gray scales indicate 
the ^^C0(2-l) and ^^C0(2-l) line emission, as indicated in each box. Contour levels increase from 0.2 in steps of 0.2 Jy/beam for 
CH3CN, from 0.24 in steps of 0.24 Jy/beam for ^^CO SMA, from 0.6 in steps of 1 Jy/beam for ^^CO SMA+30m, from 0.12 in steps 
of 0.175 Jy/beam for ^^CO SMA, and from 0.15 in steps of 0.5 Jy/beam for ^^CO SMA+30m. The cross in the bottom left of each 
box denotes the angular and spectral resolutions. 



From the continuum flux (see Table [2]), we derive a mass of 
1700 Mq for the HMC, assuming a dust temperature of 100 K 
and a dust absorption coefficient at 1.3 mm of 0.005 cm^ g"^ 
(Kramer et al. 120031) . The temperature of 100 K is a lower limit, 
because the real temperature must be higher than or equal to 
the maximum brightness temperature of the CH3CN lines mea- 
sured in the synthesized beam (-95 K). The latter corresponds 
to the surface temperature of the HMC, because the ground state 
CH3CN lines are optically thick and, as discussed by BEL05, 
the temperature inside the HMC is likely increasing. The mass 
obtained is aff'ected by a large uncertainty mostly because of 
the poorly known dust properties: for example, BEL04 adopted 
an absorption coefficient of 0.02 cm^ g"\ which would imply a 
HMC mass of only 420 M©. 

In Fig. \T3\ we show the same position- velocity plot of the 
CH3CN and ^^CO lines as in Fig. [121 this time overlaying the 
pattern outlining the region inside which emission is expected 
for a Keplerian rotating and free-falling disk. The latter does not 
take into account the line- width nor the spectral and angular res- 
olutions and has been obtained assuming that the gas velocity is 
the vector sum of a tangential component due to Keplerian ro- 
tation about a central mass, M, plus a radial component due to 



free-fall onto the same mass. Under these assumptions the ve- 
locity component along the line of sight can be expressed as 



V = ^/GM— -F ^2GM-j 
R-2 R2 



(1) 



where x and z are the coordinates, respectively, along the disk 
plane and the line of sight, and R = V^^Tz^ is the distance 
from the center of the disk. We also assume that R lies between 
the disk radius, Rq, and a minimum inner radius, Ri. The first 
term on the right hand side of the equation is the component 
due to Keplerian rotation, the second that due to free fall. The 
dashed pattern in Fig.[T2]has been obtained by plotting the max- 
imum and minimum velocities V for z varying across the disk, 
i.e. from - ^JRl - x^ and -\- y^^-^, taking into account that 
the region R < Ri is forbidden. In our case, a satisfactory fit is 
obtained for M = 330 M©, Ro = ^ (or 0.11 pc), and Ri = Of! 15 
(or 0.0055 pc). Note that the signature of (pseudo-)Keplerian 
rotation is the "butterfly" shape of the plot, determined by the 
two "spurs" of emission at about ±3'' and ±4 km s"^ relative 
to the systemic velocity (-96.5 kms"^), plus the presence of 



14 



R. Cesaroni et al.: Dissecting a hot molecular core: The case of G31.41+0.31 



110 



100 



E 



> 




2 0-2-4 

^NE-sw (arcsec) 

Fig. 13. Overlay of the ^^C0(2-l) (only SMA data; contours) 
and CH3CN(12-11) K=4 position- velocity plots along the di- 
rection with P.A.=68° passing through the HMC. The offset is 
measured from the phase center, positive toward NE. Contour 
levels increase from 0.2 in steps of 0.2 Jy/beam for CH3CN and 
from 0.1 in steps of 0.1 Jy/beam for ^^CO. The vertical and hori- 
zontal solid lines denote the position of the continuum peak and 
the systemic velocity, respectively. The thick dashed pattern en- 
compasses the region inside which emission is expected from 
a disk with radius of 3'' (or 0.11 pc) undergoing free-fall and 
Keplerian rotation about 330 M©. The cross in the bottom left 
denotes the angular and spectral resolutions. 



high- velocity emission at zero offseO- It is worth stressing that 
the contribution of the ^^C0(2-l) line is crucial to outline such 
a pattern, which cannot be recognized from the sole CH3CN 
emission. This explains why BEL05 managed to fit their CH3CN 
data assuming a flat or solid-body rotation curve, and Girart et 
al. (12009b , from a number of CH3OH lines, found that the rota- 
tion velocity increases with radius. 

One may wonder whether a HMC as massive as several 
100 Mq may be undergoing Keplerian rotation. As suggested by 
Cesaroni et al. (12007b and demonstrated by Beltran et al. (201 1 ), 
these massive, large rotating cores are to be considered transient 
toroidal structures feeding a cluster of YSOs rather than stable 
circumstellar accretion disks. In fact, the latter are stabilized by 
the central star(s), whose mass is greater than that of the disk, 
whereas the former are dynamically dominated by the gas mass 
and hence are short-lived. Although this scenario may be true 
in most cases, G31.41 might represent an exception. We specu- 
late that this HMC could contain a large number of stars tightly 
packed in the central region. In this case the stars could have a 
stabilizing eff'ect analogous to that of a single point-like object 
located at the HMC center, similar to the previously mentioned 
model by Berlin & Lodato ([T999b . 

With all this in mind, the pseudo-Keplerian pattern recog- 
nized in Fig. [13] suggests that the HMC mass should be com- 
parable to the total mass of the embedded stars that are tightly 



^ The peak in the contour plot at C and ~84 km s~^ is likely due to 
an unidentified line of a core tracer, whereas the broad wing of emis- 
sion extending up to ~70 km s~^ is genuine high-velocity ^^C0(2-l) 
emission. 



packed at the center. Is this scenario plausible? Indeed, in the 
case of GlO.62-0.38, a massive star forming region with a lumi- 
nosity of 9.2 X 10^ Lo, Sollins et al. (2005 ) suggest that within 
a radius of 0.03 pc, several O stars with a total mass of 175 M© 
have formed at the center of a flattened disk. In our case, the ex- 
istence of (at least) two high-mass YSOs close to the HMC cen- 
ter and separated (in projection) by only 0^.^19 or 1480 AU has 
been proved by Cesaroni et al. (2010), who detected the free- 
free continuum emission. In terms of mass, these two YSOs rep- 
resent only the tip of the iceberg, because the total mass of the 
associated cluster is dominated by (undetected) low-mass stars. 
Assuming e.g. a Miller & Scalo (1979) mass function, one can 
estimate the cluster mass from that of the most massive star. For 
this we can use the value of 20-25 M© computed by Osorio et 
al. C2009') from their model fit and obtain a total mass of the clus- 
ter of ~ 10^ Mq. The same result is obtained by fixing the total 
luminosity of the cluster to the value of 2 x 10^ L© quoted by 
Osorio et al. (.2009 ). In all likelihood ~ 10^ M© is an upper limit 
to the real stellar mass, because the cluster mass function may be 
highly incomplete in such a small region. However, this estimate 
indicates that the dynamical mass of -330 M©, obtained from the 
fit in Fig. \T3\ is probably dominated by the stellar mass, which 
would lend support to the application of Berlin & Lodato (.1999D 
model to the case of G3 1 .41 . 

Finally, we note that the toroid interpretation is compatible 
with the hourglass-shaped morphology of the magnetic field (see 
Girart et al. 2009), whose symmetry axis is directed SE-NW, 
because the latter coincides with the rotation axis of the toroid. 
However, a caveat is in order. Girart et al. (120091) find a signifi- 
cant correlation between rotation velocity and radius in the core, 
using diff'erent tracers. In particular, they conclude that in the 
region sampled by their observations (lying between ~0'.'5 and 
IV 6), velocity increases with distance from the HMC center (see 
their Fig. 4). This result seems inconsistent with our hypothesis 
that the core is undergoing pseudo-Keplerian rotation, because 
in this case the rotation velocity should decrease with radius. 
Perhaps this discrepancy can be explained by the presence of 
the magnetic field. Higher angular resolution observations are 
needed to sample the velocity field inside the core and thus ob- 
tain a reliable, direct measurement of the rotation curve. 

4.3. The CH3CN "8-shaped" structure 

The last aspect we will discuss in the context of the out- 
flow/toroid controversy, is the double-peaked, "8-shaped" struc- 
ture observed in the CH3CN K=2 map at the systemic velocity 
(see Fig.O. The same morphology is seen in all components up 
to K=6, i.e. up to excitation energies of -300 K, whereas the 
^=8 map (energy of 513 K) presents a single peak at the HMC 
center (Fig. [3]). In Sect. 13.11 we interpreted these facts in terms 
of opacity and temperature gradients. Indeed, BEL05 found that 
the CH3CN rotational temperature peaks toward the HMC center 
and from the ratio between our CH3CN and CH3^^CN (12-11) 
K=2 data we calculate an optical depth at the peak of the line in 
the range 8-77 across the HMC. The interplay between opacity 
and temperature can explain the existence of a dip at the HMC 
center, but does not justify the lack of circular symmetry in the 
low-^ maps. Here, we wish to find an explanation for the "8- 
shaped" feature and check whether this can better fit into the 
outflow or toroid model. 

As a basis for our discussion, in Fig. [14] we present an over- 
lay of the blue- and red-shifted CH3CN K=4 emission on the 
bulk emission in the same line. In this way, one is comparing the 
high- with the low- velocity emission. 
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Fig. 14. Maps of the blue- (solid contours) and red-shifted 
(dashed) emission in the CH3CN(12-11) K=4 line overlayed 
on the bulk emission (gray scale) in the same line. The 
integration was performed over the velocity ranges 88.9- 
91.3 km s-^ 100.9-103.3 km s-^ and 95.2-96.7 km s-^ respec- 
tively. Contour levels increase from 0. 15 in steps of 0.27 Jy/beam 
for the blue- and red-shifted gas, and from 0.2 in steps of 
0.2 Jy/beam for the bulk emission. The synthesized beam of the 
images is drawn in the bottom left. The two white points mark 
the positions of the two compact free-free sources detected by 
Cesaroni et al. (l2010i) . 



In the outflow scenario, a naive interpretation of this figure 
is that the high- velocity gas is leaking from the HMC through 
the axis of a "donut-like" structure seen edge-on, corresponding 
to the "8-shaped" feature in the map. In this case, the two peaks 
would coincide with the maxima of column density across the 
"donut". Albeit plausible, this interpretation has a problem. The 
large CH3CN opacity should prevent the detection of two dis- 
tinct peaks, because the line brightness is independent of column 
density. Therefore, instead of the 8-shaped feature, one should 
see an elongated structure perpendicular to the bipolar outflow 
and peaking at the center of it (i.e. at the HMC center). 

Explaining the presence of the 8-shaped structure appears to 
be a problem also in the toroid scenario. In fact, if the toroid 
is seen edge-on, the emission at the systemic velocity should 
peak in between the blue- and red-shifted emission. However, 
if the toroid is inclined with respect to the line of sight, the 
nearest and farthest sides of it are seen displaced from one an- 
other on the plane of the sky, symmetrically displaced with re- 
spect to the center. This is indeed what we see in Fig. [141 where 
the two peaks lie along the (projection of the) rotation axis and 
are equally ofl'set from the center. The situation depicted in this 
figure resembles that of the ring undergoing Keplerian rotation 
about the low-mass system GG Tau (see Fig. 6a of Guilloteau et 
al. 11999 ). Clearly, this comparison is inappropriate, because the 
scales of the two objects differ by more than an order of mag- 
nitude. Moreover, G31.41 is probably more similar to a "puffy 
pancake" rather than a "donut" or ring, because gas and dust 
are present at all radii, as demonstrated by the continuum and 
high-energy line emission peaking at the center of the HMC (see 
Figs. [2] and O. Note, however, that the emission from the lower 



energy lines may be arising mostly from the colder outer region 
and hence be confined in a ring. Finally, G3 1 .41 , unlike GG Tau, 
is embedded in a dense pc-scale envelope, which complicates 
the line and continuum radiative transfer. Despite all these dif- 
ferences, the comparison between the two objects is intriguing 
and may qualitatively explain the features observed in G3 1.41. 

Further support for this interpretation is obtained by com- 
paring the case of G31.41 to that of a source closer in mass than 
GG Tau. This is the Herbig Ae star MWC758, which is sur- 
rounded by a circumstellar disk, as demonstrated by the model 
fit of Isella et al. (12010b . A comparison between the channel 
maps presented by these authors (see their Fig. 2) with those 
in Fig. [4] reveals a surprising similarity: in both cases the shape 
of the emission is circular and peaks at opposite positions with 
respect to the center in the blue- and red- shifted channels, but 
turns into an 8-shaped structure close to the systemic velocity. 
This fact may be the fingerprint of a rotating disk in G3 1.41, as 
wen as in MWC 758. 

5. Summary and conclusions 

With the present study we have performed a deep analysis of the 
geometrical structure and kinematics of the HMC in the high- 
mass star- forming region G31.41. Thanks to new high angular 
resolution SMA images and complementary IRAM 30-m maps, 
we have investigated both the HMC, through the CH3CN(12-1 1) 
transition, and the lower density surroundings, through ^^CO and 
^^CO. We could thus search for a possible outflow powered by 
the star(s) embedded in the HMC. In particular, our main goal 
was to shed light on the nature of the NE-SW velocity gradient, 
previously detected in the HMC by various authors, and possibly 
distinguish between rotation and expansion. 

Our new data confirm the presence of the velocity gradi- 
ent basically in all HMC tracers observed, and indicate that a 
red-blue symmetry along the same axis is seen also in the CO 
isotpomers. We were unable to detect any bipolar outflow along 
the SE-NW direction, in contrast with the findings of Olmi et 
al. (1996). The latter was probably an artifact caused by very 
limited uv sampling and lack of zero- spacing information. 

We find a hint of infall in the HMC, consistent with the re- 
sults of Girart et al. (2009). A comparison between the CH3CN 
and CO maps suggests that these species are tracing the same 
phenomenon. Especially the ^^C0(2-l) line emission appears to 
be the obvious prosecution of the CH3CN(12-1 1) emission on a 
larger scale. 

We have discussed whether the velocity gradient observed 
in G31.41 is a compact bipolar outflow or a rotating, geomet- 
rically thick, toroidal structure. The former hypothesis implies 
outflow parameters typical of high-mass stars with luminosities 
in excess of 10^ L©, whereas previous estim ates by Osorio et 
al. (^2009 ) and Cesaroni et al. C 1994al l20T0b indicate a signifi- 
cantly lower luminosity for G31.41. Moreover, such an outflow 
would be much more compact than typical bipolar outflows de- 
tected in high-mass star- forming regions and the dynamical time 
scale would be an order of magnitude shorter than that needed to 
form HMC spacies such as CH3CN. It is also worth noting that a 
core undergoing expansion seems difficult to reconcile with the 
detection of infall. 

The composite CH3CN and ^^CO position- velocity plot can 
be explained if the gas is undergoing both (pseudo-)Keplerian 
rotation and infall, suggesting that this massive object could be 
dynamically stabilized by a compact cluster of YSOs tightly 
packed inside a few 1000 AU from the center. In this scenario 
the ^^CO emission could be tracing the outer regions, where the 
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gravitational field is dominated by the gas and the rotation veloc- 
ity tends to a constant value. Albeit intriguing, this interpretation 
remains speculative and could be proved only by producing syn- 
thetic CH3CN and CO maps and position- velocity plots to be 
compared with the data. This numerical eff'ort goes beyond the 
purposes of the present study. 

We conclude that the case of G31.41 could be a scaled-up 
version of lower-mass YSOs, such as GG Tau and MWC758, 
where circumstellar disks have been detected and successfully 
modeled. We stress that whatever the interpretation, the obser- 
vations of this object have produced an important finding for the 
study of high-mass star formation. We obtained iron-clad evi- 
dence for the existence of a velocity gradient across the HMC 
and demonstrated that this gradient is not due to multiple cores 
unresolved in the beam and moving at diff'erent velocities, but 
to the gas undergoing a smooth drift in (line-of- sight) velocity 
from the one end to the other of the core. If, as we believe, 
one is dealing with a rotating toroid, then this indicates that 
star formation in this HMC may proceed in a similar way to 
low-mass stars, namely through circumstellar, centrifugally sup- 
ported disks "hidden" inside the densest interiors of the toroid 
itself. If, instead, the outflow interpretation were correct, such 
a sharp, neat velocity gradient would strongly suggest the ex- 
istence of a very efl'ective focusing mechanism for the outflow 
and the best candidate as a collimating agent would be a (yet 
undetected) circumstellar disk. Therefore, regardless of whether 
the G31.41 HMC is undergoing expansion or rotation, the role 
of disks in the mechanism of high-mass star formation seems to 
be strengthened by our findings. The advent of ALMA will be 
crucial to detect deeply embedded disks in HMCs. 
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